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ABSTRACT 

In  this  thesis,  a  new  scheme  of  encoding  orthogonal  MFSK  signals  with  convolutional  codes 
is  proposed.  Instead  of  separating  the  code  and  the  signals  as  is  done  conventionally,  here  the 
m  outputs  of  the  convolutional  encoder  of  rate  k/m  are  used  to  select  the  M  =  2m  orthogonal 
FSK  signals.  This  combined  coded  system  is  called  trellis-coded  modulation  MFSK 
(TCM/MFSK). 

The  free  distance  of  TCM/MFSK  codes  is  found  computationally.  Asymptotic  coding  gains 
(ACG)  of  several  TCM/MFSK  systems  relative  to  uncoded  MFSK  systems  having  the  same 
bandwidth  efficiency  (bits/s/Hz)  are  calculated. 

Both  analytic  descriptions  and  natural  mapping  implementations  of  TCM/4-FSK  with  rate  1/2 
and  TCM/8-FSK  with  rate  2/3  for  several  constraint  lengths  of  the  convolutional  codes  are 
given.  The  analytic  description  of  TCM/MFSK  is  also  obtained  computationally. 
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I.    INTRODUCTION 

The  model  of  the  digital  communication  system  under  consideration  is  shown  in 
Figure  1.1.  The  input  to  the  convolutional  encoder,  {aj,  a^  =  0,  1  is  a  sequence  of  k 
independent  bits.  The  channel  encoding  operation  is  combined  with  M-ary  Frequency 
Shift  Keying  (MFSK)  modulation  as  in  Trellis  Coded  Modulation  (TCM)  schemes.  A 
convolutional  encoder  is  used  to  generate  the  underlying  code.  MFSK  orthogonal 
signaling  is  used  to  achieve  an  acceptable  performance  with  a  minimum  power 
requirement. 

The  set  of  orthogonal  signals  used  is  denoted  as: 


{S(t)},  x=l,2, M 


(1.1) 


These  signals  are  generated  by  a  frequency  synthesizer  which  is  driven  by  a  set  of  real 
numbers 


n(t) 


0(t> 


Convolution*-! 


r-fc/B 


^(t) 


r(t) 


Dmxlalmtor 


^(t) 


•i.k 


Figure  1 . 1   System  Model 


{*},  jc  =  1,2, M  (12) 

These  real  numbers  have  a  specific  relationship  with  the  incoming  bit  stream  fa}  and  the 
memory  bits  of  the  encoder.  For  an  encoder  that  has  2V  states,  the  constraint  length  of 
the  convolutional  code  is  defined  as  v.  The  constraint  length  v  is  equal  to  the  memory 
of  the  encoder  (shift  register  stages). 

The  M   =   2m  distinct  symbols  are  represented  by  M  orthogonal  sinusoidal 
waveforms  as: 

Sx(t)=Acos[2ir(fc  +  (x-\)Af)t],0<t<ZT  (1-3) 

where 

x  =  1,  2,  M 

T:  signaling  interval 

fc:  carrier  frequency 

Af:  the  width  of  the  frequency  slot. 
The  orthogonal  signal  set  is  characterized  by  equal  signal  energy  E  given  by: 

T 

E=ls2x(t)dt  =  ^I  (1-4) 


I 


The  orthogonality  condition  is  satisfied  when  the  cross-correlation  coefficient  obeys  the 
relationship 


T 

prUt)Sjw=  c/J*  (i.5) 


E  ifi=j 


MFSK  is  generated  by  subdividing  a  frequency  interval  into  M  distinct  frequency  slots. 
Each  slot  has  a  width  Af.  The  minimum  frequency  separation  between  adjacent  signals 
is  obtained  from  (1.5)  [Ref.  1] 

A/    =1  (1.6) 

->  mm        t 

Thus,  the  minimum  one-sided  bandwidth  occupied  by  these  M  orthogonal  signals  is: 

B  =  ™  (1.7) 

T 

It  is  assumed  that  the  transmission  medium  introduces  zero  mean,  additive  white 

N 
Gaussian  noise  (AWGN)  n(t)  with  power  spectral  density  _£ .  The  received  signal  is 

2 
corrupted  by  AWGN  as: 

rx(t)=sx(t)+n(t)  (1.8) 

This  signal,  rx(t),  is  demodulated  by  a  bank  of  either  matched  filters  or  correlators. 
When  the  initial  phase  can  be  estimated  by  the  receiver,  the  demodulation  is  coherent, 
otherwise  it  is  noncoherent. 

The  analog  outputs  (ux(t)},  x  =  1,  2,  ...  M  of  the  demodulator  are  then  sent  to  the 
Viterbi  decoder  which  performs  a  soft  decision  decoding  of  the  MFSK  signal.  The 
sequence  {aj  at  the  output  of  the  Viterbi  decoder  is  the  maximum-likelihood  estimate  of 


the  input  sequence  {aj.  In  this  thesis,  only  coherent  TCM/MFSK  systems  are 
considered. 

The  remaining  chapters  in  this  theisis  are  organized  as  follows. In  Chapter  II ,  an 
overview  of  convolutions!  coding  and  TCM/MFSK  are  presented.  The  results  from  the 
computer  program  written  to  find  the  free  distance  of  TCM/MFSK  are  given.  Asymptotic 
coding  gain  (ACG)  is  defined  and  a  list  of  ACGs  is  also  given.  Chapter  HI  deals  with  the 
design  of  the  encoder  for  TCM/MFSK.  Both  the  analytic  expressions  and  the  natural 
mapping  rules  for  certain  schemes  are  examined.  Although  all  the  signals  in  the  M-ary 
orthogonal  constellation  have  equal  distances  from  each  other,  it  is  shown  that  the 
transmitter  complexity  depends  on  the  signal  assignment.  Basic  rules  are  given  for  the 
minimal  transmitter  complexity.  In  Chapter  IV,  several  encoders  are  designed  on  a  step 
by  step  basis.  TCM/MFSK  is  further  discussed  and  some  conclusions  are  drawn  in 
Chapter  V. 

The  computer  program  to  find  the  free  euclidean  distance  of  TCM/MFSK  is  given 
in  Appendix  A.  Another  program  to  find  the  analytic  expression  of  TCM/MFSK  is  given 
in  Appendix  B. 


H.   OVERALL  REVIEW 

In  1948,  Shannon  showed  that  the  main  limitation  to  transmission  over  a  noisy 
channel  is  set  by  a  quantity  called  channel  capacity.  Capacity  of  a  noisy  channel  is  the 
largest  rate  at  which  information  can  be  transmitted  reliably.  In  other  words,  if  the  data 
source  rate  is  less  than  the  channel  capacity,  proper  encoding  and  decoding  techniques 
enable  us  to  communicate  over  a  noisy  channel  with  any  arbitrary  error  rate.  Otherwise, 
reliable  communication  is  not  possible. 

Codes  are  divided  into  two  categories  according  to  their  function.  Codes  for  source 
encoding  attempt  to  reduce  the  data  source  rate  at  the  price  of  tolerable  distortion.  Codes 
for  channel  encoding  turn  a  noisy,  constrained  channel  into  a  reliable,  unconstrained 
channel.  The  application  of  channel  encoding  is  called  Forward  Error  Correcting  Coding 
(FEC).  Most  communication  systems  use  FEC  coding  systems.  FEC  coding  is  generally 
achieved  by  redundancy  adding  and  noise  averaging.  Codes  are  classified  according  to 
their  structure  as  block  codes  and  convolutional  codes. 

In  this  thesis,  channel  encoding  and  MFSK  modulation  operations  are  combined 
together  to  take  advantage  of  soft  decision  decoding.  Convolutional  coding  is  used  for 
this  purpose.  It  may  be  possible  to  use  another  code  as  the  outer  code  to  improve 
performance  further. 

The  basic  distinction  between  convolutional  codes  and  block  codes  is  that 
convolutional  codes  do  not  have  a  fixed  block  length.  These  codes  are  generated  by 


passing  the  input  sequence  through  a  system  of  linear  shift  registers  and  modulo-2 
adders.  A  convolutional  code  is  characterized  by  a  ratio  called  the  code  rate.  The  code 
rate  is  the  ratio  of  the  number  of  input  lines  to  the  number  of  output  lines  of  the  encoder. 
Thus,  the  code  rate  is 

r  =  A  (2.1) 

m 

where  k  is  the  number  of  inputs  and  m  is  the  number  of  outputs  of  the  encoder.  The 
performance  of  the  encoder  depends  strongly  on  the  number  of  states  2V  where  the 
constraint  length  v  is  equal  to  the  memory  of  the  encoder. 

The  basic  utility  of  convolutional  codes  comes  from  the  fact  that  these  codes  can 
be  decoded  by  the  Viterbi  decoding  algorithm.  The  Viterbi  algorithm  is  an  efficient  way 
of  sequence  demodulating  in  the  presence  of  noise  [Ref.  2].  The  Viterbi  decoding 
algorithm  is  equivalent  to  maximum-likelihood  decoding  and  thus  optimum  for  equally 
likely  messages  [Ref.  3]. 

When  the  encoding  and  modulation  are  done  separately,  the  basic  parameter 
affecting  the  code  performance  is  the  Hamming  distance,  which  is  the  number  of  places 
in  which  any  two  sequences  differ.  If  the  encoding  and  the  modulation  are  combined  into 
one  step,  the  Euclidean  distance  governs  the  performance  of  the  coded  scheme.  The 
Euclidean  distance  between  two  sequences  c  and  d  is  defined  as 

where  cl  and  dt  are  the  elements  of  c  and  d,  respectively. 


At  the  receiver  end,  if  the  Viterbi  decoder  deals  with  the  quantized  demodulator 
outputs  which  are  the  best  estimates  of  the  transmitted  signal  (i.e.,  the  number  of 
quantization  levels  is  2),  it  is  called  hard-decision  decoding.  If  the  Viterbi  decoder  works 
on  the  sequence  of  information  on  which  a  final  decision  can  be  made  about  the 
transmitted  signals  (i.e.,  the  number  of  quantization  levels  is  more  than  2),  the  method 
is  called  soft-decision  decoding.  Soft-decision  decoding  has  more  coding  gain  than 
hard-decision  decoding  [Ref.  1]. 

Trellis  coded  modulation  (TCM)  proposed  by  Ungerboeck  combines  encoding  and 
modulation  to  design  a  code  for  a  particular  geometry  of  the  signal  constellation.  Thus, 
the  Euclidean  distance  becomes  the  main  parameter  and  soft-decision  decoding  is  always 
used.  The  trellis  encoder  is  composed  of  two  parts.  The  first  part  is  the  convolutional 
encoder  with  v  bits  of  memory.  This  determines  the  number  of  the  states  2V  of  the 
encoder.  The  second  part  is  the  mapper.  This  part  is  memory  less.  Coded  digits  select  the 
signal  from  the  constellation  and  move  the  encoder  to  the  next  state.  Uncoded  digits  only 
select  signals  causing  parallel  transitions.  The  mapper  maps  the  incoming  coded  and 
uncoded  digits  to  the  M  signals.  The  mapping  rule  is  defined  by  a  function  which  must 
be  nonlinear  in  order  to  achieve  a  coding  gain  [Ref.  4]. 

The  asymptotic  coding  gain  is  the  ratio  of  EJN0  required  for  coded  operation 
relative  to  that  required  for  uncoded  operation  with  the  same  bandwidth  efficiency  in  the 
limit  of  large  EJN0.  Although  actual  coding  gain  is  a  more  realistic  indicator  of  system 
performance,  it  is  significantly  more  difficult  to  obtain  than  is  ACG,  and  ACG  does  give 


an  indication  of  the  performance  that  can  be  expected  for  AWGN  channels.  The  ACG 
of  TCM/MFSK  is  defined  as: 


ACG  =  lO\og(-dlj2E)        for  TCM/4-FSK  (2.3a) 
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ACG  =  101og(IrfjL/2£)        for  TCM/S-FSK  (2.3b) 

where 

dfree!  free  Euclidean  distance  of  the  coded  scheme  (the  smallest  of  the  Euclidean 
distances  between  any  two  coded  paths) 

E     :  average  energy  of  the  coded  signal 

Decoding  is  done  using  the  Viterbi  algorithm.  It  is  basically  a  search  of  the  most 
likely  path  through  the  trellis.  Sometimes  the  chosen  path  may  not  be  the  transmitted  path 
due  to  the  effect  of  noise.  Then  an  error  event  is  said  to  have  taken  place.  An  error  event 
is  the  part  of  the  path  which  diverges  from  the  correct  path  and  remerges  at  a  later  time. 
The  free  distance  of  a  TCM  scheme  is  the  minimum  Euclidean  distance  of  all  error 
events.  Thus,  by  definition,  free  distance  plays  an  important  role  in  determining  the 
performance  of  any  TCM  scheme. 

A  general  analytic  solution  for  the  free  distance  can  not  be  obtained;  hence,  it  is 
found  computationally.  A  computer  program  is  given  in  Appendix  A  to  find  the  free 
distance  of  any  TCM/MFSK.  When  the  signal  set  is  orthogonal,  the  distance  between  any 
pair  of  M  signals  is: 


dtj  =  yJlE         for  all  ij  (2-4) 

This  is  also  the  minimum  distance  of  any  orthogonal  signal  constellation.  In  the  computer 
program  this  distance  is  normalized  to  unity  regardless  of  the  number  of  the  signals  in 
the  constellation.  Thus, 

dtj  =  sf2E    =  1         for  all  i,j  (2-5) 

The  number  of  states  in  the  trellis  diagram  is  2V.  The  number  of  the  transitions  that 
diverge  from  a  common  state  or  remerge  into  a  same  state  is  2k.  Since  the  code  rate  is 
r  =  k/m  with  m  _>_  k  +  l,M  =  2m  signals  are  assigned  to  2k+v  transitions.  The  figure 
of  merit  to  compare  different  coding  schemes  is  ACG  given  by  (2.3).  The  reference 
uncoded  constellation  is  the  one  which  has  the  same  bandwidth  efficiency  as  the  coded 
scheme.  Convolutional  coding  does  not  effect  the  spectrum  of  the  transmitted  signal 
[Ref.  5].  Hence,  the  bandwidth  efficiency  R/B  may  be  taken  as  a  reference  for 
comparison  purposes.  Here  R  is  the  information  rate  in  bits  per  second  and  B  is  the 
minimum  one-sided  bandwidth  in  Hz  given  by  (1.7).  For  an  uncoded  M-ary  orthogonal 
signaling,  where  M  =  2k'  the  information  rate  is 


R  =  *!  (2.6) 

T 


where  k'  =  log2M  and  T  is  the  symbol  duration.  Also  B  =  M/T,  hence: 


—  =  —  =  _  (2  7) 

B       M       2« 

When  a  convolutional  encoder  with  rate  r  =  k/m  is  used,  the  information  rate  is  R  = 
k/T,  but  the  bandwidth  is  expanded  to  B  =  2m/T.  Thus,  the  coded  bandwidth  efficiency 
is 

*  =  JL  (2.8) 

B       2m 

ACG  is  calculated  for  the  coded  schemes  for  which  (2.7)  and  (2.8)  are  equal.  For 
example,  a  rate  xh  TCM/4-FSK  with  k  =  1 ,  m  =  2  has  a  bandwidth  efficiency  of  V* 
b/s/Hz.  Uncoded  16-FSK  with  k'  =  4  has  the  same  bandwidth  efficiency. 

A  search  with  the  program  given  in  Appendix  A  shows  that  the  squared  free 
distance  of  the  best  codes  is 

4f  =L  +  1  (2.9) 

IE 

where  L  =  |_v/kj  ,  and  [_•  J  denotes  the  integer  part.  The  results  are  tabulated  in  Table 

2.1.  Note  that  the  bandwidth  efficiency  for  all  TCM/4-FSK  and  TCM/8-FSK  schemes 

is  14  b/s/Hz. 

The  program  input  consists  of  a  2V  x  2k  matrix  N.  The  rows  of  N  are  the  transitions 

from  the  states  of  the  trellis  respectively.  An  example  of  a  trellis  diagram  and  the 

corresponding  matrix  N  are  given  in  Figure  2.1.  The  states  are  labeled  as  1,  2,  ,T. 

The  first  row  of  N  consists  of  the  labels  of  the  destinations  of  the  transitions  from  the 

first  state.  The  second  row  for  the  transitions  from  the  second  state,  so  on. 
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Figure  2.1    A  Trellis  Diagram  and  its  Matrix 
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TABLE  2.1  ASYMPTOTIC  CODING  GAINS  OF  TCM/4-FSK  AND  TCM/8-FSK 
WITH  RESPECT  TO  UNCODED  16-FSK 


Code  Rate 
(r) 

Modulation 
(MFSK) 

Constraint 
Length  (v) 

d]rtt 

IE 

ACG  (dB) 

1/2 

4-ary 

2 

3 

-1.25 

1/2 

4-ary 

3 

4 

0.0 

1/2 

4-ary 

4 

5 

0.97 

1/2 

4-ary 

5 

6 

1.76 

1/2 

4-ary 

6 

7 

2.43 

1/2 

4-ary 

7 

8 

3.01 

1/2 

4-ary 

8 

9 

3.52 

1/2 

4-ary 

9 

10 

3.98 

1/2 

4-ary 

10 

11 

4.39 

2/3 

8-ary 

2 

2 

0.0 

2/3 

8-ary 

3 

2 

0.0 

2/3 

8-ary 

4 

3 

1.76 

2/3 

8-ary 

5 

3 

1.76 

2/3 

8-ary 

6 

4 

3.01 

2/3 

8-ary 

7 

4 

3.01 

2/3 

8-ary 

8 

5 

3.98 

2/3 

8-ary 

9 

5 

3.98 

2/3 

8-ary 

10 

6 

4.77 
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m.    ENCODER  DESIGN 

The  function  of  the  encoder  is  to  map  the  input  binary  stream  into  M  orthogonal 
signals.  This  may  be  viewed  as  a  sliding  window  method  of  encoding  a  binary  data 
stream  into  a  sequence  of  M  real  numbers.  If  the  number  of  input  bits  into  the  encoder 
is  k  and  the  code  constraint  length  is  v,  then  the  channel  input  Sx(t)  will  depend  on  n  = 
k  +  v  bits.  The  v  bits  in  memory  determine  one  of  2V  states  of  the  encoder  while  2k 
transitions  are  associated  with  k  input  bits. 

The  MFSK  waveform  is  given  by  (1.3)  and  it  is  the  output  of  a  voltage  controlled 
oscillator  (VCO)  driven  by  a  set  of  real  numbers  (voltage  levels).  Hence,  the  encoding 
operation  may  be  implemented  by  mapping  n  bits  into  the  set  of  real  numbers  {x},  x  =  1, 
2,  ...  M.  Because  x  is  a  real-valued  function  of  n  bits,  it  may  be  written  as  a  sum  of 
products  of  input  bits  ^  where  a^  €{0,1}  or  b;  where  b;  €{-1,1}  and  b;  =  l-2a; 
[Ref.  6].  Thus, 


n 
n  V 

z=l         v. 

J>i  (3.1) 

n 

+  i]%= 1 4"MA +  dn.J\  b2 K 

h>j>i 


where 

x  (bi,  b2,  bj  G  {x},  x  =  1,  2,  M  and  d,,  d^,  d^, 
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are  the  coefficients  that  we  want  to  find,  and  bi5  br  \\  E{-1,1}. 
Following  the  discussion  in  [Ref.  4],  we  note  that  there  are 


n\    n\    n 


0  +  l+2 


*G  -  1" 


(3.2) 


coefficients.  Let  x  denote  a  column  vector  with  2n  elements  which  are  all  the  values  that 
x  can  assume.  Also,  let  d  denote  a  column  vector  whose  elements  are  the  unknown 
coefficients.  Let  B  be  a  2n  x  2n  matrix  where  each  row  represents  the  2°  values  taken  by 
all  the  products  of  b;'s  in  (3.1).  Thus,  (3.1)  may  be  written  in  matrix  form  as  follows: 


x=Bd 


where 


x  = 


41,1,1, 1) 

4-1,1,1, 1) 


4-1,-1,-1,....-!)     v 


2"xl 


and 

Bk  =  [1  D!  b2b12 b12 J 

is  the  k*  row  of  2n  x  2n  matrix  B.  Also, 

cF  =  [do  dj  d2  d12 d12 J 


(3.3) 
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where  the  superscript  T  denotes  matrix  transpose.  We  note  that  B  is  a  2n  x  2n  square 
matrix  and  B  =  BT,  hence  B  is  symmetric.  In  other  words,  it  is  an  orthogonal  matrix 
[Ref.  4].  Therefore: 

BBT  =  2nI  (3-4) 

where  I  is  2°  x  2n  identity  matrix.  Thus, 

Bl  =  —BT  (3.5) 


2" 


and  the  solution  of  (3.3)  is 


d  =  ±BTx  (3.6) 

2" 

A  computer  program  listing  is  given  in  Appendix  B  to  calculate  the  coefficients  in 
d.  The  program  requires  the  number  of  inputs  k,  the  constraint  length  v,  and  the  vector 
jc.  The  elements  of  x  are  from  the  set  {x},  x  =  l,  2,  ...  M.  Thus,  x  is  the  signal 
assignment  of  the  corresponding  trellis  diagram. 

Once  the  relation  between  the  input  stream  {bj}  and  {x},  x=l,  2,  ...  M  is  known, 
the  encoder  is  realized  according  to  logic  variables  fa}.  In  other  words,  a  natural 
TCM/MFSK  mapper  may  be  implemented  if  the  analytic  description  of  the  encoder  is 
available.  In  an  orthogonal  MFSK  signal  constellation,  all  signals  have  equal  euclidean 
distance  from  each  other  as  stated  in  (2.3).  Thus,  the  distance  is  of  no  importance  when 
the  signals  in  x  are  assigned.  Instead,  care  should  be  given  to  avoiding  a  catastrophic 
code  (a  limited  number  of  errors  causes  an  infinite  number  of  errors  at  the  Viterbi 
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decoder  output)  and  intersymbol  inference  should  be  kept  in  mind.  Another  important 
concern  may  be  the  complexity  of  the  encoder.  If  the  three  rules  below  are  followed, 
generally  one  ends  up  with  a  good  code: 

1 .  All  elements  of  x  should  be  equiprobable. 

2.  The  first  2m  elements  must  consist  of  the  elements  ofx  =  l,  2,  M. 

3.  The  frequency  differences  must  be  kept  constant  giving  the  maximum  signal 
difference  to  bx. 

These  rules  are  derived  from  a  computer  search  of  many  signal  assignments.  Rule  1 

ensures  to  keep  the  optimality  of  decoding  which  is  maximum-likelihood  and  optimal  for 

the  equally  likely  messages.  All  three  rules  are  also  intended  to  avoid  catastrophic 

conditions  and  reduce  Doppler  interference.  If  we  end  up  with  a  consistent  mapping  rule, 

this  ensures  that  TCM/MFSK  is  not  catastrophic.  Rule  3,  with  the  other  two  rules, 

guarantees  that  adjacent  transitions  have  the  maximum  frequency  separation  which 

reduces  Doppler  interference. 

A  frequency  difference  corresponds  to  the  absolute  difference  between  any  two 

elements  of  {x},  x  =  l,  2,  ...  M  whose  associated  trellis  transitions  differ  only  one  bit. 

For  example,  the  frequency  difference  due  to  bj  is 

«,  =  \x{bvb2i ^ bkJ-x{b{b2 -bj ,bkJ\  (3.7) 

When  all5j'sj  =  l,2,...nare  kept  constant,  d  has  a  minimum  number  of  nonzero  elements 
and  transmitter  complexity  is  minimal  [Ref.6]. 
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A  list  of  the  codes  for  TCM/4-FSK  and  TCM/8-FSK  with  several  constraint  lengths 
is  given  at  the  end  of  this  chapter.  In  signal  assignment  for  4-ary  codes,  the  frequency 
differences  are  kept  constant  in  a  periodic  fashion,  as  follows: 

51=2,  52  =  1,  53  =  2,  54  =  1  (3.8) 

On  the  other  hand,  for  8-ary  codes,  the  following  frequency  differences  are  used. 

{5„  /  =  1,2,...10}  =  {4,2,1,4,2,1,4,2,2,4}  (3.9) 

For  TCM/8-FSK,  k  is  2  and  the  codes  with  constraint  length  v  up  to  8  are  considered. 
Hence,  the  maximum  n  =  2  +  8  =  10.  Because  the  channel  waveform  depends  on 
n  =  k  +  v  bits,  TCM/8-FSK  waveforms  with  v  =  8  will  be  a  function  of  10  bits.  If  the 
constraint  length  v  is  taken  to  be  2,  TCM/8-FSK  waveforms  are  determined  by 
n  =  k+  v  =  4  bits.  Thus,  only  8U  62,  53  and  54  will  be  considered. 

Once  we  have  an  analytic  expression  for  the  code  we  may  find  the  natural  mapping 
rule  in  terms  of  the  logic  variables  a,,  a^  E{0,1}.  We  note  that  ai  =  VHl-bj)  and  that 
multiplication  with  b(  variables  corresponds  to  an  exclusive  or  function  in  logic  variables 
ai.  The  natural  mapping  is  the  process  of  mapping  the  output  of  the  convolutional  encoder 
to  one  of  the  M  signals  according  to  a  given  rule.  For  example,  Gray  encoding  is  a  kind 
of  mapping  with  the  rule  that  adjacent  signals  have  a  change  of  only  one  bit.  The  natural 
mapping  is  defined  in  logical  terms,  not  in  physical  quantities;  hence,  the  constants  in 
analytic  expression  have  no  meaning  in  the  natural  mapping. 

The  codes  in  logic  variables  corresponding  to  analytic  expressions  in  Table  3.3  are 
listed  in  Table  3.4.  The  variables  {Zj},  j  =  1,  2,  ...  m  are  the  inputs  to  the  natural 
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frequency  mapper.  (Obviously,  Zj's  can  also  be  viewed  as  the  outputs  of  the  convolution^ 
encoder.)  The  mapping  rules  for  TCM/4-FSK  and  TCM/8-FSK  are  found  by  examining 
the  signal  assignment  of  each  scheme  and  are  given  in  Table  3.1  and  Table  3.2, 
respectively. 

TABLE  3.1   TCM/4-FSK  NATURAL  MAPPING  RULE 


Zi 

z2 

Frequency  Slot 

0 

0 

1 

0 

1 

2 

1 

0 

3 

1 

1 

4 

The  analytic  expression  and  consequently  the  natural  mapping  rule  is  a  function  of 
signal  assignment.  Thus  it  is  possible  to  find  other  codes  by  changing  the  signal 
assignment  within  the  rules  given  previously.  In  other  words,  the  code  for  a  given 
constraint  length,  code  rate,  and  modulation  is  not  unique. 
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TABLE  3.2   TCM/8-FSK  NATURAL  MAPPING  RULE 


zt 

^ 

Z, 

Frequency  Slot 

0 

0 

0 

1 

0 

0 

1 

2 

0 

1 

0 

3 

0 

1 

1 

4 

1 

0 

0 

5 

1 

0 

1 

6 

1 

1 

0 

7 

1 

1 

1 

8 

TABLE  3.3   ANALYTIC  DESCRIPTIONS  OF  THE  CODES 

Code 
Rate  (r) 

TCM/MFSK 

Constraint 
Length  (v) 

Analytic 
Description 

1/2 

4-ary 

2 

x=2.5-0.5b2-b1b3 

1/2 

4-ary 

3 

x=2.5-0.5b2b4-b,b3 

1/2 

4-ary 

4 

x=2.5-0.5b2b4-b1b3b5 

1/2 

4-ary 

5 

x=2.5-0.5b2b4b6-b1b3b5 

1/2 

4-ary 

6 

x=2.5-0.5b2b4b6-b1b3b5b7 

1/2 

4-ary 

7 

x=2.5-0.5b2b4b6b8-b1b3b5b7 

2/3 

8-ary 

2 

x=4.5-0.5b3-b2-2b,b4 

2/3 

8-ary 

3 

x=4.5-0.5b3-b2b5-2b1b4 

2/3 

8-ary 

4 

x=4.5-0.5b3b6-b2b5-2b1b4 

2/3 

8-ary 

5 

x=4.5-0.5b3b6-b2b5-2b1b4b7 

2/3 

8-ary 

6 

x=4.5-0.5b3b6-b2b5b8-2b1b4b7 

2/3 

8-ary 

7 

x=4.5-0.5b3b6b9-b2b5b8-2b1b4b7 

2/3 

8-ary 

8 

x=4.5-0.5b3b6b9-b2b5b8-2b1b4b7b10 
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TABLE  3.4   CODES  IN  LOGIC  VARIABLES 


Code  Rate 
(r) 

TCM/MFSK 

Constraint 
Length 

z, 

z2 

z, 

1/2 

4-ary 

2 

a,  0a-, 

*2 

- 

1/2 

4-ary 

3 

3,0a, 

a20a4 

- 

1/2 

4-ary 

4 

a,  ©a,0a< 

a-0a4 

- 

1/2 

4-ary 

5 

ai0a3©a5 

a20a40a6 

- 

1/2 

4-ary 

6 

a|0a30a50a 

7 

a20a40a6 

- 

1/2 

4-ary 

7 

a,0a30a50a 

7 

a20a40a60a 

8 

- 

2/3 

8-ary 

2 

a,0a4 

*2 

h 

2/3 

8-ary 

3 

a,0a4 

a20a5 

*3 

2/3 

8-ary 

4 

a,0a4 

a20a5 

a30a6 

2/3 

8-ary 

5 

a,0a40a7 

a20a5 

a30a6 

2/3 

8-ary 

6 

a,0a40a7 

a~0as0a>. 

^0^ 

2/3 

8-ary 

7 

a,0a40a7 

a20a50ag 

a30a60a9 

2/3 

8-ary 

8 

a|0a40a70a 

10 

a,0as0a^ 

a30a60a9 

©  denotes  exclusive  or 
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IV.    ENCODER  DESIGN  EXAMPLES 

A.      ENCODER  DESIGN  1 

In  this  example,  a  rate  V2  TCM/4-FSK  encoder  is  considered.  The  constraint  length 
v  is  chosen  to  be  3.  There  are  v  =  3  memory  elements  and  2V  =  8  trellis  states,  each 
having  2k  =  2  transitions.  Normalized  squared  free  Euclidean  distance  is  found  by  using 
the  program  in  Appendix  A,  which  is  4.  The  input  to  this  program  is  the  matrix  N.  The 
trellis  structure,  signal  assignment,  matrix  N,  and  input/state  variables  are  shown  in 
Figure  4.1. 


1.3    • 


N= 


1  2 

3  4 

5  e 

7  a 

1  2 

3  4 

5  6 

7  § 

Matrix  W 


**  I 1  **  I 1  b>    I \b* 


Trellis  /  Signal   Assignment  Input  /  State   Variables 

(The  elements  of  x) 


Figure  4. 1 
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The  rows  of  N  are  the  labels  of  the  destination  states,  respectively.  For  example, 
the  first  state  has  transitions  to  state  1  and  state  2;  hence,  the  first  row  of  N  has  the 
labels  1  and  2.  The  second  state  has  transitions  to  state  3  and  state  4,  and  the  elements 
of  the  second  row  are  3  and  4  and  so  on. 

The  signal  assignment  is  made  by  keeping  the  frequency  differences  constant; 
hence, 

5j  =2,  52  =  1,  53  =  2,  54  =  1 

For  example,  the  transition  from  state  1  to  state  1  occurs  when  the  first  frequency  is 
transmitted.  If  the  third  frequency  is  transmitted  then  the  new  encoder  state  is  the  second 
one.  The  definition  of  frequency  difference  is  stated  in  (3.7).  The  adjacent  transitions  are 
different  in  the  b,  bit  only,  thus  the  frequency  difference  d1  due  to  bj  is  2.  The  signal 
assignment  in  Figure  4.1  is  put  in  vector  form  as: 

x  =  [132431422413423  1] 
The  input  of  the  program  given  in  Appendix  B  is  the  x  vector.  The  program  also  requires 
the  number  of  inputs  k  and  the  constraint  length  v.  The  outputs  of  this  program  are  the 
coefficients  in  d  which  are  found  to  be: 

d0  =  2.5 

d13  =  -1 

d24  =  -0.5 
Thus,  from  (3.3)  the  analytic  description  of  the  code  is 
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x  =  2.5-b{b3-0.5b2b4 


0.5 


Figure  4.2   Analytic  Description  Transmitter  for  TCM/4-FSK,  8  States 

The  analytic  description  of  the  transmitter  may  be  implemented  as  in  Figure  4.2.  The 
natural  mapper  implementation  in  terms  of  the  logic  variables  ^  is  as  follows.  If  we  make 
the  change  of  variables 

x  =  2.5-yl-0.5y2 

where  v7  =  b}bs  and  v2  =  bp4,  then  since  the  multiplication  in  bj  corresponds  to  an 
exclusive-or  function  in  logic  variables,  we  can  find  the  new  variables  as 


z2  =  a2©a4 
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These  new  variables  zlf  z2  are  input  to  the  natural  mapper.  The  mapping  rule  is  found 
by  examining  the  signal  assignment.  Thus,  the  natural  mapper  implementation  may  be 
as  in  Figure  4.3  with  the  following  mapping  rule. 

TABLE  4.1 


Zl 

z2 

Frequency 

0 

0 

1 

0 

1 

2 

1 

0 

3 

1 

1 

4 

a. 


4  + 


Figure  4.3   Natural  Mapper  for  r  =  Vfc,  v  =  3,  TCM/4-FSK 
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The  bandwidth  efficiency  of  rate  Vi  TCM/4-FSK  scheme  can  be  found  from  (2.7) 


as 


JL  =  A  blslHz 
2m       4 

Uncoded  16-FSK  has  the  same  bandwidth  efficiency.  The  ACG  is  0  dB  with  respect  to 
uncoded  16-FSK  (Table  2.1). 

B.      ENCODER  DESIGN  2 

A  rate  2/3  TCM/8-FSK  with  constraint  length  v  =  4  is  introduced  in  this  example. 
There  are  4  memory  bits  and  2V  =  16  trellis  states.  The  number  of  transitions  at  each 
state  is  2k  =  4.  The  squared  free  euclidean  distance  is  found  to  be  3.  The  trellis  diagram, 
signal  assignment,  matrix  N,  and  input/state  variables  are  shown  in  Figure  4.4. 

The  bandwidth  efficiency  of  the  rate  2/3  TCM/8-FSK  encoder  is  (2.7) 

—  =  I  blslHz 
2m       4 

Uncoded  16-FSK  has  the  same  bandwidth  efficiency.  The  ACG  with  respect  to  uncoded 
16-FSK  is  1.76  dB  (Table  2.1). 
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1,5,3,7 
2,f,4,$ 
5,1,7,3 
6,2,8,4 
3,7,1,5 
4,$,2,f 

7,  3,  5, 1 
8,4,4,2 
2,6,4,8 

1,5,3,7 
6,2,8,4 
5,1,7,3 
4,8,2,6 
3,7,1,5 
8,4,6,2 
7,3,5,1 


N 


3 

4 

7 

6 

11 

12 

15 

If 

3 

4 

7 

8 

11 

12 

13 

If 

3 

4 

7 

8 

11 

12 

15 

If 

3 

4 

7 

8 

11 

12 

15 

If 

1      2 
5      6 
9      10 
13   14 
1      2 
5      6 
9     10 
13   14 
1      2 
5      6 
9      10 
13   14 
1      2 
5      6 
9      10 
13   14 


TielliM  tutrix 


i   r 

2*         »     w  ■ 


b    . .  bA  , .    b. 

2  *  * 

_»     T  —    T   — . 


Input/State   Vaxlablaa 


Tre 11  la /Signal  AMllgnnent 


Figure  4.4 
When  the  signal  assignment  in  Figure  4.4  is  put  in  a  vector  x  as 


x  =  [1537264851736284 
3715482673518462 
2648153762845173 
4826371584627351] 


and  input  to  the  program  in  Appendix  B,  the  d  coefficients  are  found  to  be 
do  =  4.5 
d14  =  -2 
d25  =  -1 
d36  =  -0.5 
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Hence,  the  analytical  description  of  the  set  {x},  x  =  1,  2,  8  in  terms  of  b;  is 

x  =  A.5-2bxbA-bj}s-Q.5b^6 

The  above  equation  may  be  written  as 

x  =  4.5-2v1-y2-0.5v3 

where 

yi  =  bjb4 

y2  =  bp5 
y3  =  bp6 

Therefore,  z,  variables  in  terms  of  at  are 

z7  =  a1®a4 

z2  =  a2®a5 

z3  =  a3®a6 
with  the  following  mapping  rule  (Table  4.2).  The  natural  mapper  implementation  of  rate 
2/3  TCM/8-FSK  with  16  states  is  shown  in  Figure  4.5. 
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TABLE  4.2 


Z3 

z2 

Zl 

Frequency 

0 

0 

0 

1 

1 

0 

0 

2 

0 

1 

0 

3 

1 

1 

0 

4 

0 

0 

1 

5 

1 

0 

1 

6 

0 

1 

1 

7 

1 

1 

1 

8 

Figure  4.5   Natural  Frequency  Mapper  for  r  =  2/3,  v  =  4,  TCM/8-FSK 
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C.      ENCODER  DESIGN  3 

A  rate  Vi  TCM/4-FSK  encoder  with  32  states  is  considered.  The  constraint  length 
is  v  =  log2  32  =  5.  Each  state  has  2  transitions.  The  squared  free  Euclidean  distance  is 
6;  ACG  with  respect  to  uncoded  16-FSK  is  1.7  dB  (Table  2.1).  The  matrix  N  which  is 
the  definition  of  the  trellis  diagram  is  given  in  (4.1)  (following  page). 

The  signal  assignment  keeps  the  following  frequency  differences  constant 

5j=2,  52  =  1,  53=2,  54  =  1,  65  =  2,  56  =  1 

The  vector  x  is 

[1324142  2  413423131421324 
jc=    423124132413423113243142 
4231241331421324] 

The  above  signal  assignment  yields  the  following  analytic  expression 

x  =  2.5-b2b,b6-bxb,b5 

and  the  z,  variables  are 

Zj  =  a}®a3®a5 

z2  =  a2®a4®a6 
with  the  mapping  rule  given  in  Table  3.1.  The  natural  mapper  implementation  is  shown 
in  Figure  4.6. 
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AT  = 


1  2 
3  4 
5  6 
7  8 
9  10 
11  12 
13  14 
15  16 
17  18 
19  20 
21  22 
23  24 
25  26 
27  28 
29  30 
31  32 
1  2 
3  4 
5  6 
7  8 
9  10 
11  12 
13  14 
15  16 
17  18 


19  20 

21  22 

23  24 

25  26 

27  28 

29  30 

_  31  32 

(4.1) 
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z, 


Frequency 
Mapper 


S(t) 


Figure  4.6  Natural  Frequency  Mapper  rate  =  Vz,  v  =  5,  TCM/4-FSK 


D.      ENCODER  DESIGN  4 

A  TCM/8-FSK  encoder  with  32  states  is  given  to  make  the  subject  clearer.  The 
rate  is  r  =  2/3;  hence  k  =  2,  m  =  3.  The  constraint  length  v  is  5  and  ACG  is  1.76  dB. 
The  matrix  N  is  given  in  (4.2). 
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N  = 


12  3  4 
5  6  7  8 
9  10  11  12 
13  14  15  16 
17  18  19  20 
21  22  23  24 
25  26  27  28 
29  30  31  32 
12  3  4 
5  6  7  8 
9  10  11  12 
13  14  15  16 
17  18  19  20 
21  22  23  24 
25  26  27  28 
29  30  31  32 
12  3  4 
5  6  7  8 
9  10  11  12 
13  14  15  16 
17  18  19  20 
21  22  23  24 
25  26  27  28 
29  30  31  32 


(4.2) 
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In  this  frequency  assignment  the  following  differences  are  kept  constant  (3.9) 

5,  =4,  52  =  2,  52  =  1,  54  =  4,  55  =  2,  86  =  1,  57=4 

These  constraints  yield  the  following  x  vector. 

[15372648517362843715 
48267351846226481537 
62845173482637158462 
x =  73515173628415372648 
73518462371548266284 
51732648153784627351 
4  8  2  6  3  7  15] 

The  x  vector  has  2v+k  =  128  elements.  The  following  analytic  expression  is  obtained  by 
the  signal  assignment  above 

x  =  4.5-0.5^6-Z?2Z75-2^4Z?7 

This  analytic  expression  corresponds  the  following  zt  variables  with  the  rule  given  in 
Table  3.2 

Zj  =  aj®a4®a7 

z2  =  a2®a5 

z3  =  a3@a6 
Finally,  the  natural  mapper  can  be  implemented  as  shown  in  Figure  4.7. 
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Figure  4.7  Natural  Frequency  Mapper,  r  =  2/3,  v  =  5,  TCM/8-FSK 
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V.   DISCUSSION  AND  CONCLUSIONS 

The  objective  of  this  thesis  was  to  investigate  trellis  coding  of  M-ary  orthogonal 
signals.  The  power  efficiency  of  this  class  of  signaling  is  attractive  to  many 
communication  systems  when  the  environment  is  power-limited.  An  example  of  such  an 
environment  is  the  generation  of  low  probability  of  intercept  (LPI)  signals.  MFSK  is  used 
as  an  example  of  orthogonal  signaling,  but  the  discussion  holds  for  other  orthogonal 
signalings  such  as  Pulse  Position  Modulation  (PPM). 

A  method  of  encoding  orthogonal  signals  by  using  a  smaller  constellation  with 
convolutional  coding  for  a  given  bandwidth  efficiency  was  considered.  The  free 
Euclidean  distance  of  TCM/MFSK  codes  plays  an  important  role  when 
maximum-likelihood  Viterbi  soft-decision  decoding  is  used  on  an  AWGN  channel.  A 
computer  program  has  been  written  to  fmd  the  free  Euclidean  distance.  A  computer 
search  of  rate  Vi  TCM/4-FSK  and  rate  2/3  TCM/8-FSK  shows  that  the  normalized 
squared  free  Euclidean  distance  increases  linearly  with  the  constraint  length,  as; 


dfree  .    V 


=  m+i 


IE        "Jfc" 

where  v  is  the  constraint  length,  k  is  the  number  of  input  bits,  and  |_*J  is  the  "floor" 
operation  (the  integer  part).  The  normalized  distance  between  any  pair  of  signal  points 
is  taken  to  be  unity;  i.e., 
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dtj  =  yJlE    =  1         for  all  ij    (i  *f) 
The  investigation  shows  that  an  arbitrarily  small  error  rate  can  be  achieved  with  a 
TCM/MFSK  scheme  by  using  a  larger  v. 

ACG  comparisons  are  made  with  respect  to  uncoded  16-FSK  having  the  same 
bandwidth  efficiency.  ACGs  of  4.39  dB  and  4.77  dB  for  TCM/4-FSK  and  TCM/8-FSK 
with  constraint  length  10  are  observed,  respectively.  This  means  TCM/MFSK  can 
maintain  better  performance  while  using  the  same  bandwidth  than  uncoded  MFSK. 

Convolutional  coding  with  rates  k/m  are  proposed  to  be  used  as  the  underlying 
code.  This  choice  brings  several  advantages.  First,  efficient  soft-decision  Viterbi 
decoders  can  be  used.  Second,  since  rate  k/m  convolutional  codes  are  used,  the  signal 
space  is  expanded  as  much  as  2m/2k  in  terms  of  convolutional  encoder  input-output 
relations.  It  is  shown  that  more  coding  gain  is  achievable  with  this  method  than  the 
conventional  approach.  Doubling  the  signal  space  provides  the  greatest  coding  gain  while 
requiring  a  minimum  increase  in  system  complexity  [Ref.  7].  Finally,  an  outer 
coding  structure  may  be  used  for  further  improvements  in  system  performance. 

A  computer  program  was  written  to  find  the  analytic  description  of  the  codes.  Once 
the  analytic  description  is  known,  the  natural  mapper  implementation  is  straightforward. 
The  basic  rules  of  frequency  assignment  for  minimal  transmitter  complexity  were 
discussed.  The  key  is  to  keep  the  frequency  differences  constant.  It  is  possible  to  find 
different  codes  by  changing  the  signal  assignment  only.  Encoder  design  examples  are 
worked  out  which  are  also  significant  as  another  way  of  convolutional  code  design. 
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Although  ACG  gives  an  idea  about  the  performance  on  an  AWGN  channel,  the 
actual  coding  gain  of  trellis  coded  orthogonal  signaling  on  both  AWGN  channels  and 
fading  channels  for  both  coherent  and  noncoherent  demodulation  is  required  to  complete 
our  understanding  of  the  performance  that  can  be  expected  from  this  type  of  system 
under  realistic  conditions. 
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APPENDIX  A 


%%%%%%%%%%%%%     PROGRAM  DISTANCE     %%%%%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 
% 

%  This  program  is  written  in  Matlab  to  find  the  free  distance  of  a 
%  trellis.  The  transitions  from  each  state  are  saved  in  a  matrix  N. 
%  This  matrix  is  saved  in  a   filename. m  file. 

%  The  first  row  of  the  matrix  N  consists  of  the  #'s  of  the  states  which 
%  can  be  reached  from  state  1  in  1  trans it ion. And  the  second  row  is  for 
%  the  second  state... (The  first  state  is  l,not  0) .The  program  takes  care 
%  of  the  rest. The  parallel  transtions  are  not  considered. The  program  asks 
%  the  name  of  the  input  file. The  output  of  the  program  is  the  squarred 
%  free  euclidean  distance  of  the  trellis  described  by  the  matrix  N 
% 
% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear; 

%  Input  initilaziation 

r»input (' Enter  the  input  file  name      ','s'); 

N=eval(r); 

%  s  :  #  of  the  states 

%  t  :  #  of  the  transitions  at  each  state 
[s, t] -size(N) ; 

%  Initilazation  of  the  state  transition  matrix 
B=zeros (s, s) ; 
for  i=l:s 

for  j-l:t 

z-N(i,j); 
B (i,  z) -z; 
end 
end 
for  i=l:s 

for  j-l:s 

if  j«B(i,j) 
B(i,j)-1; 
else 

B(i,j)-100; 
end; 
end; 
end; 

%  Calculation  of  the  shortest  distances 
for  k»l:s 

for  i-l:s 

if  B(i,k)-  =  100 
for  j=l:s 

length=B(i,k) +B(k, j) ; 
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if  length<B(i, j) 

B{i, j) -length; 
end; 
end; 
end; 
end; 
end 

%  Calculation  of  the  free  distance  w.r.t.  all  zero  codeword 
for  j=l:s 

if  B(l,j)~l 

DF(l,j)-B(j,l)+B(l,j); 
DF(1,1)«100; 
else 

DF(l,j)-100; 
end 
end 

%  Squarred  free  distance 
d2free=min(DF) 

%  Pairwise  Comparision  for  Checking 

%  If  the  result  of  the  pairwise  comparision  does  not 

%  agree  the  trellis  setup  is  wrong 

dfree=100; 

yeni=B; 

for  i=l:s 

if  B(i,i)--1 
yeni (i,i) =100; 
for  j-l:s 

free=yeni (i, j) +yeni ( j , i) ; 
if  free<dfree 
dfree=free; 
end 
end 
end 
end 
for  i=l:s 

for  j=l:s 

if  yeni (i, j) --1 
yeni (i,j) =100; 
for  x»l:s 

if  yeni (i,x) »»1 
for  k-l:s 

if  B(j,k)«B(x,k) 
free»B(j  ,k)  +1; 
if  free<dfree 
dfree-free; 
end 
end 
end 
end 
end 
end 
end 
end 
dfree 
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APPENDIX  B 


%%%%%%%%    PROGRAM  CODE   %%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% 

%  This  program  finds  the  coefficients  for  the  analytic  desciption  of  a 
%  TCM/MFSK  scheme. The  inputs  to  the  program  are  the  #  of  the  inputs   k, 
%  the  constraint  length  v,  and  a  file.  The  input  file  contains  a  vector  X 
%  which  is  the  signal  assignment  of  the  trellis. The  number  of  the 
%  elements  is  equal  to  the  total  number  of  the  transitions. 
% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
clear 

%  Input  initilazion 

m=input ('Enter  the  number  of  input  bits        ','s') 
r=input ('Enter  the  costraint  length  ','s') 

t -input (' Enter  the  name  of  the  input  file     ','s') 

%  k  :  #  of  the  input  bits 

%  v  :  the  constraint  length 

%  X  :  the  signal  vector 

%  n  :  #  of  input  bits  plus  the  number  of  the  memory  bits 

k=eval (m) ; 
v=eval (r) ; 
X=eval(t) ; 

n=k+v; 

%  Creating  Hadamard  matrix 

B=hadamard(n) ; 

%  Calculating  coefficients 

D=B'*X'/2An; 

%  Determining  the  place  of  the  coeffients 

d-find(D') 
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